120° Helical Magnetic Order in the Distorted Triangular Antiferromagnet a-CaCr 2 04 
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a-CaCr2C>4 is a distorted triangular antiferromagnet. The magnetic Cr 3+ ions which have spin- 
3/2 and interact with their nearest neighbors via Heisenberg direct exchange interactions, develop 
long-range magnetic order below Tjv = 42.6 K. Powder and single-crystal neutron diffraction reveal 
a helical magnetic structure with ordering wavevector k = (0, ~ 1/3,0) and angles close to 120° 
between neighboring spins. Spherical neutron polarimetry unambiguously proves that the spins lie 
in the ac plane perpendicular to k. The magnetic structure is therefore that expected for an ideal 
triangular antiferromagnet where all nearest neighbor interactions are equal, in spite of the fact 
that Q-CaCr2 04 is distorted with two inequivalent Cr 3+ ions and four different nearest neighbor 
interactions. By simulating the magnetic order as a function of these four interactions it is found 
that the special pattern of interactions in a-CaCr2 04 stabilizes 120° helical order for a large range 
of exchange interactions. 



PACS numbers: 75.25.-fz, 61.05.F-, 75.30.Et, 75.40.Cx 
I. INTRODUCTION 

Frustrated magnets, characterized by competing mag- 
netic interactions, continue to generate much research 
interest. In these systems it is not possible to satisfy 
all magnetic interactions simultaneously, as a result the 
ground state can be highly degenerate leading to exotic 
physical states e.g. spin liquid behaviour and chiral or- 
der. The simplest frustrated system is the triangular 
lattice antiferromagnet where all magnetic interactions 
between nearest neighbors ( J nn ) are equal. In their pio- 
neering work, Anderson and Fazekas suggested that the 
ground sta te i s a spin liquid with no long-range mag- 
netic order However, recent theoretical work implies 
that at T = it develops long-range order^ with a helical 
structure, where the spin moments on nearest neighbors 
point 120° with respect to each other. Interactions be- 
tween the planes can stabilize this ground state at finite 
temperatures. 

Among real triangular lattice materials the 120° 
structure has been found in CsFc(SC>4)2, RbFc(S04)2, 4 
VX 2 (X=C1, Br, I) El CuCrCbP Ag 2 Ni0 2 pl and 
RbFe(MoC>4) 2 P Unfortunately, these compounds only 
exist in polycrystalline form or as small single-crystals 
limiting the possibilities of experimental investigation. 
Additional terms in the spin Hamiltonian can favour dif- 
ferent magnetic structures or even destroy long range or- 
der. For example single-ion anisotropy favors collinear 
order (CuFeC>2p cv-NaMnC^ 1 -); frustrated interlayer in- 
teractions suppress order (AgCr0 2 p^ NaCr02^) and 
strong exchange striction drive the system a way from 
120° helical magnetism (a-NaMn0 2 P CuCrSa ^ 1 ) . 



Departures from ideal triangular crystal symmetry 
typically lead to spatially anisotropic exchange interac- 
tions sometimes accompanied by orbital ordering and 
resulting in the frustration being lifted. In most cases 
the distortion reduces the dimensionality, so that one of 
the three nearest neighbor exchange interactions (J nn i) is 
stronger than the other two (J nn 2), and together the in- 
traplanar interactions produce antiferromagnetic chains 
with frustrated interchain interactions. The magnetic 
structure is often helical but where the angle between 
nearest neighbours along the chain takes a value between 
120° and 180°. 

Examples of such compounds are CS2CUCI4 where 
"Air^/^nni = 0.33 and the angle between intrachain neigh- 
bours is 170° while for a-NaMn0 2 the presence of sin- 
gle ion anisotropy along with the ratio Jni^/^nn i = -44 
gives rise to collinear antiferromagnetic orderP^ESl A 
third example is CuCr02, a multiferroic compound where 
a small lattice distortion accompanies the transition to 
long-range magnetic order. This along with substantial 
next nearest neighbour and frustrated interplane interac- 
tions, res ults in angles between nearest neighbors of 118° 
and 120°™ 

In this paper we will investigate a-CaCr 2 04 which 
is a triangular lattice antiferromagnet belonging to the 
delafossite arisotypeP^ Although this compound is also 
distorted from ideal triangular symmetry (orthorhom- 
bic space group Pmmn), the distortion does not lower 
the dimensionality and is of a different type to that ex- 
plored previously. Figure [T] shows the crystal structure. 
There are two inequivalent Cr 3+ ions (spin-3/2) which 
together form distorted triangular layers in the be plane, 
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that long range magnetic order develops below TV =43 
K with ordering wavevector k = (0, 0.3317(2), 0) and 
spins lying in the ac plane. The magnetic structure is 
indirectly inferred to be helical rather than sinusoidal 
so as to constrain the size of chromium moment within 
physical limits, and the angles between nearest neighbor 
spins are almost 120°. Using symmetry arguments it is 
shown that this compound is not ferroelectric although 
quadratic magneto-electric effects are possible. 

We have grown the first single-crystals of a-CaCr 2 C>4 
and in this paper we investigate the magnetic proper- 
ties using heat capacity, DC magnetic susceptibility, neu- 
tron and X-ray powder diffraction, neutron single-crystal 
diffraction and spherical neutron polarimetry. We show 
that despite the fact that the crystal structure is distorted 
from ideal triangular symmetry, the magnetic structure 
is that expected for an undistorted triangular lattice anti- 
ferromagnet. The magnetic ordering wavevector is close 
to the commensurate value k = (0, 1/3,0) and the spins 
form a helical structure in the ac plane with ~ 120° be- 
tween nearest neighbors. By simulating the magnetic 
structure as a function of the four independent nearest 
neighbor exchange interactions we show how the specific 
type of distortion in a-CaCr 2 04 is able to give rise to the 
observed highly symmetric magnetic order. 



II. EXPERIMENTAL DETAILS 



FIG. 1. Crystal structure and magnetic interactions of a- 
CaCr2 04. (a) Crystal structure showing the triangular layers 
formed from two inequivalent Cr06 octahedra colored blue for 
Cr(l)06 and red for for Cr(2)06- The layers are separated 
by Ca 2+ ions, (b) A single triangular layer showing only the 
Cr 3+ ions. The four inequivalent nearest-neighbour intralayer 
Cr 3+ -Cr 3+ distances determined from the E9 powder data at 
2.1 K (see Table [i]): red d chl = 2.882 A; blue d ch2 = 2.932 A; 
green d zz i = 2.914 A; purple d ZZ 2 = 2.910 A. The correspond- 
ing exchange interactions are J c hi, Jch2, Jzzi and J ZZ 2, (see 
Discussion section). The dashed lines indicate the boundaries 
of the structural unit cell. 



that are stacked along a. The Cr 3+ ions are in an oc- 
tahedral environment and the three d electrons fill the 
three t 2g orbitals resulting in quenched orbital moment 
and negligible anisotropy. Neighboring CrOe octahedra 
are edge-sharing leading to direct overlap of the i 2g or- 
bitals and this along with the short Cr 3+ -Cr 3+ distances 
(~ 3 A) result in direct antiferromagnetic exchange in- 
teractions. Superexchange via oxygen is suppressed be- 
cause the Cr-O-Cr angles are ~ 90° . There are a total of 
four different nearest neighbor distances giving four in- 
equivalent exchange interactions which together form the 
complex pattern shown in Fig. [ljb) . Interplane distances 
are much larger (~ 5.5 A) and only weak superexchange 
interactions are expected between the layers. Recently 
the magnetic properties of cv-CaCr 2 04 were investigated 
by Chapon et alW^ Using powder samples they show 



Both polycrystalline and single-crystal samples of a- 
CaCr 2 04 were synthesized. The polycrystalline samples 
were produced by a solid state reaction. A 1:1 molar ratio 
of CaC03 and Cr 2 C>3 was mixed thoroughly, pressed into 
a pellet and annealed in N 2 atmosphere at 1100°C for 12 
hours. The powder was then quenched in liquid nitrogen 
to prevent the formation of impurity phases. The single- 
crystal growth was performed using a high-temperature 
optical floating zone furnace at the Crystal Laboratory, 
Helmholtz Zentrum Berlin fur Materialien und Energie 
(HZB), Berlin, Germany, details of the technique will be 
reported elsewherePH The single-crystal used for neutron 
scattering experiments is plate-like with shiny flat sur- 
faces perpendicular to the a-axis, it has a weight of 340 
mg and a size of 3 x 7 x 8 mm 3 . This crystal consists of 
three crystallographic twins rotated by 60° with respect 
to each other about the shared a-axis. Smaller single- 
crystal pieces were used for the bulk properties measure- 
ments. 

Bulk properties measurements were performed at the 
Laboratory for Magnetic Measurements, HZB. Heat ca- 
pacity was measured on a small single-crystal sample of 
15.6 mg using a Physical Properties Measurement System 
(PPMS), Quantum Design over the temperature range of 
2.1 < T < 300 K. High temperature DC magnetic sus- 
ceptibility (300 < T < 1000 K) was measured in a 1 
T field on an unoriented crystal using the PPMS with 
vibrating sample magnetometer and oven options. Low 
temperature (2.1 < T < 300 K) susceptibility at 1 T 
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and magnetisation up to 5 T were measured both paral- 
lel and perpendicular to the a-axis (x a || and XaJ_) using a 
superconducting quantum interference device (Magnetic 
Property Measurement System, Quantum Design). The 
single-crystal had a weight of 1.81 mg and was fixed with 
GE varnish to a plastic sample stick. Both field cooled 
and zero field cooled susceptibility were measured but no 
significant difference was observed. 

X-ray powder diffraction was measured on the ID31 
high-resolution powder-diffraction beamline at the Euro- 
pean Synchrotron Radiation Facility (ESRF), Grenoble, 
France. A wavelength of A = 0.39983 A (E = 31.18 
keV) was used, and data were collected at a number 
of temperatures in the range 10 < T < 260 K. Neu- 
tron powder patterns were measured on two instruments. 
First on the D20 high-intensity two-axis diffractometer at 
the Institute Laue Langevin (ILL), Grenoble, France. A 
Ge monochromator was used to select a wavelength of 
A = 1.869 A, and data were collected over the temper- 
ature range of 10 < T < 42.6 K in 1 K steps. Mea- 
surements were also performed on the E9 fine-resolution 
powder diffractometer at the BER II reactor, HZB. A 
Ge monochromator selected a wavelength of A = 1.799 A 
and data were collected at 2.1 K and 50 K. 

Single-crystal neutron diffraction was performed on the 
large (340 mg) crystal using the E5 four-circle single- 
crystal diffractometer at the BER II reactor, HZB. A 
pyrolytic graphite (PG) monochromator selected an in- 
cident wavelength of A = 2.36 A. The instrument is 
equipped with a two-dimensional position-sensitive 3 He 
detector consisting of 36 x 36 pixels and for each Bragg 
reflection a 36-step scan of the crystal angle (w) was per- 
formed. A total of 616 magnetic and 423 nuclear Bragg 
reflections were collected at 6 K for a total counting time 
of five minutes each. 

Spherical neutron polarimetry was performed on the 
same single-crystal using the TASP triple axis spectrome- 
ter at the Paul Scherrer Institute (PSI), Switzerland. The 
wavelength A = 3.2 A was selected by a PG monochro- 
mator and benders were used to polarize and analyse the 
neutron beam along the vertical z-axis. All measure- 
ments took place at a sample temperature of 1.5 K. The 
MuPAD option was installed to control the polarisation 
direction of the incoming neutron beam (Pi) and also to 
analyse the scattered beam along any final polarisation 
direction (Pf). The beam polarisation was measured to 
be 0.883(10) for all directions of neutron polarization. 
The conventional right-handed Cartesian coordinate sys- 
tem was used to describe the polarization with the x axis 
parallel to the wavevector transfer Q, y in the horizon- 
tal scattering plane perpendicular to Q, and z vertical. 
There are a total of 36 cross sections that can be mea- 
sured for every Bragg reflection, since both P; and Pf 
can be along ±x, ±y and ±z. For most Bragg peaks only 
18 cross sections were measured by setting the incoming 
beam polarisation to the positive value, while for a few 
peaks all 36 cross sections were measured to check chi- 
rality. Since the measured intensities were rather small, 



the background was also measured for all cross sections 
close to each Bragg peak. The polarisation matrices 
were calculated from the background corrected cross sec- 
tions: Vij = (a +i ^ +j -<T +i ^-:>)/(a +i ^ +: > +a +i ^- j ), this 
method has the advantage of canceling out positioning 
errors and any wavevector dependent attenuation since 
only the ratio of the intensities matters. A total of 23 
reflections were measured with the crystal oriented in 
three different scattering planes (h, k, 0), (h, k, 3fc/2) and 
(h,k,—3k/2). All cross sections were counted for three 
minutes each. 



III. RESULTS 

A. Bulk properties measurements 

The specific heat of a-CaCr204 as a function of tem- 
perature is shown in Fig. [2] A sharp peak is observed at 
T N = 42(1) K in agreement with Ref.HOl It is attributed 
to the onset of long-range antiferromagnetic order. No 
other sharp features were observed in this temperature 
range indicating the absence of further magnetic or struc- 
tural transitions. 
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FIG. 2. Heat capacity showing the magnetic phase transition 
at Tn = 42(1) K. The inset shows the same data plotted on 
a log-log scale. 

Figure [3^a) shows the high temperature susceptibility 
(300 < T < 1000 K), which reveals Curie-Weiss behav- 
ior. The data in the temperature range 800 K to 1000 K 
were fitted to the Curie- Weiss law: 



where 6cw is the Curie- Weiss temperature and C is the 
Curie constant measured in units of cm 3 mol _1 K _1 . xa is 
a temperature independent background, which accounts 
for the paramagnetic Van Vleck susceptibility and dia- 
magnetic core susceptibility as well the background from 
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sample stick and Al-cement (which is less than 0.5% of 
the measured susceptibility at 1000 K). The fitted value 
of the Curie-Weiss temperature is Ocw = —564(4) K 
and the effective moment which is given by \/8C/^b was 
found to be fi g = 3.68(l)/ie- This is close to the ex- 
pected spin only value 3.87/ib, assuming that the orbital 
angular momentum is completely quenched. The aver- 
age value of the exchange interactions in the triangular 
layers can be estimated from the Curie- Weiss tempera- 
ture using the mean-field approximation and is given by 
ksQ = 5(5+ l)z Jmcan/3, where z is the number of near- 
est neighbours (6 for the triangular lattice) and the in- 
terlayer interactions have been neglected. The extracted 
average intralayer coupling in a-CaCr204 is antiferro- 
magnetic with strength J moall = 6.48(5) meV. 
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FIG. 3. DC magnetic susceptibility per Cr 3+ ion measured 
on a single-crystal sample in a 1 T magnetic field, (a) left 
axis: high temperature susceptibility of unoriented sample 
(filled data points), right axis: inverse susceptibility (open 
data points), the red line is the Curie- Weiss fit. (b) Low 
temperature magnetic susceptibility; green data points, Xa,\\\ 
black points, XaX- The vertical arrow shows the inflection 
point at Tat. Right axis: inverse susceptibility (xZl) displayed 
by the open data points. 

The susceptibility at low temperatures (2.1 < T < 300 
K) for applied field parallel and perpendicular to a is 



shown in Fig. [3jb). The data deviate from the Curie- 
Weiss law and show a broad maximum at 50 K for both 
crystal orientations. This feature can be explained by the 
onset of low dimensional antifcrromagnetic correlations 
in the triangular plane. The susceptibility decreases sud- 
denly at the same temperature where heat capacity shows 
a phase transition. The large drop reveals the onset of 
long-range antiferromagnetic order and the best estimate 
for the transition temperature is Tn = 43 K. 



B. Crystal structure 

The crystal structure of cv-CaCr204 was refined from 
powder neutron diffraction data (measured on E9 and 
D20) along with X-ray powder diffraction data (mea- 
sured on ID31) collected above the antiferromagnetic 
phase transition at 50 K (Fig. [1]) . The structural refine- 
ment was performed using the FullProf softwareP All 
the peaks could be indexed in the space group Pmm n in 
agreement with previous powder measurements^^!!] anc j 
no detectable impurity was observed. There is a broad 
bump around 29 = 25° in the neutron diffraction pat- 
tern that is absent from the X-ray pattern which, due to 
the insensitivity of X-rays to magnetic moments, must 
be of magnetic origin. This provides further evidence for 
two-dimensional magnetic correlations above T/v at 50 K 
in agreement with the temperature of the maximum in 
the susceptibility (Fig. [3](b)) . The high resolution syn- 
chrotron data show that Bragg peaks with (0, k, I) are 
broadened (this was not observable in the neutron data 
because of the poorer resolution). This could be due to 
small planar rotations between adjacent layers perpen- 
dicular to the a axis. 

Table [I] lists the fitted atomic positions from the dif- 
ferent refinements. Our values are in agreement with 
those from the previous room temperature measurement 
of Ref. 19. The quality of fit is given by the i?p-factor 
values: R F =J2 \ F OBS --Fcalc|/E l-^OBsl, where F OBS 
and Fqalc are the observed and calculated structure fac- 
tors. They are: 0.035 for ID31 at 50 K, 0.037 for D20 
at 47 K, 0.035 for E9 at 50 K and 0.031 fr E9 at 2.1 K. 
The refined model confirms that a-CaCr 2 04 consists of 
distorted triangular layers with two symmetry inequiv- 
alent Cr 3+ ions in the be plane. There are a total of 
four intraplanar nearest neighbor distances varying be- 
tween 2.882 A and 2.932 A at 2.1 K, which correspond 
to four inequivalent nearest neighbour direct exchange 
paths (see Fig. |ljb) ) . The interactions follow a complex 
pattern, where equivalent exchange paths form two dif- 
ferent zig-zags and two different chains in the triangular 
plane. Cr-O-Cr angles vary between 91.0° and 95.7° 
confirming that the superexchange interactions via oxy- 
gen are weak while the large interplane distances which 
are 5.529 A (Cr(l)-Cr(l)) and 5.379 A (Cr(2)-Cr(2)) 
suggest weak interplane interactions. 

The temperature dependence of the crystallographic 
parameters were determined by sequential refinement of 
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FIG. 4. Powder diffraction patterns along with Rietveld re- 
finements. The open red symbois are the data, the black curve 
is the fit, the green bars indicate the nuclear Bragg peak po- 
sitions and the lower blue line gives the difference between fit 
and data, (a) X-ray powder diffraction measured on ID31 at 
50 K, the strongest peak (2,0,0) at 3° is excluded from refine- 
ment, (b) Neutron powder diffraction measured on D20, ILL 
at 50 K, the vertical arrow indicates the broad bump which is 
of magnetic origin and therefore missing on the X-ray pattern. 



the data collected on ID31 over a temperature range be- 
tween 10 K and 300 K. No significant change in the rel- 
ative atomic positions was observed, and in particular 
there is no sudden change in the structure at the Neel 
temperature. For comparison see the atomic positions at 
2.1 K and 50 K determined from E9 refinement in Table 
HJ Furthermore the lattice parameters were found to vary 
smoothly as shown in Fig. [5] The b and c lattice constants 
increase gradually with increasing temperature, while the 
a lattice constant is smallest at 70 K and shows negative 
thermal expansion at lower temperatures. There is no in- 
dication of any symmetry change over this temperature 
range, even a subtle distortion can be ruled out since 
the linewidth remains constant between 10 K and room 
temperature. This is in contrast to CuCr02p^ which un- 
like a-CaCr 2 04 has perfect triangular symmetry at high 
temperatures and undergoes a structural distortion at 



TABLE I. Refined crystallographic parameters using space 
group Pmrnn. The lattice parameters were determined from 
the high resolution ID31 X-ray powder data at 50K. The 
atomic position listed are; first row, ID31 powder data at 
50 K (Rf — 0.035); second row, refinement of D20 neutron 
powder data at 47 K (Rf = 0.037); third row, refinement of 
E9 powder data at 50 K (Rp — 0.035); fourth row, E9 powder 
data at 2.1 K (R F = 0.031). 



Unit cell dimensions: 



Atom Site 



a = 11.0579(2) A 
b = 5.8239(2) A 
c= 5.0553(1) A 



Ca(l) 


2a 


3/4 


3/4 


0.3512(4) 

0.3521(12) 
0.3532(12) 


Ca(2) 


2b 


1/4 


3/4 


0.0385(4) 
0.0349(10) 
0.0384(11) 
0.0410(11) 


Cr(l) 


4d 


1/2 


1/2 


1/2 


Cr(2) 


4/ 


0.4932(1) 
0.4921(4) 
0.4917(5) 
0.4931(6) 


1/4 


0.0046(4) 
0.0050(9) 
0.0059(11) 
0.0057(12) 


0(1) 


4/ 


0.4022(5) 
0.4022(4) 
0.4023(4) 
0.4024(4) 


1/4 


0.3365(14) 
0.3394(16) 
0.3397(8) 
0.3410(8) 


0(2) 


4/ 


0.5904(5) 
0.5891(3) 
0.5901(4) 
0.5906(4) 


1/4 


0.6825(14) 
0.6771(6) 
0.6758(8) 
0.6762(8) 


0(3) 


8ff 


0.5989(3) 
0.6007(2) 
0.6003(2) 
0.6003(2) 


0.4996(8) 
0.4963(5) 
0.4955(6) 
0.4966(6) 


0.1665(12) 
0.1653(5) 
0.1647(7) 
0.1657(7) 



the Neel temperature. 

Single-crystal diffraction on the E5 instrument was 
used to collect the nuclear Bragg peak intensity at 6 
K. The Racer program was used to integrate the reflec- 
tions, as described in Ref.[23l Since our single-crystal has 
a large mosaic spread (~ 2°) with several small grains 
at slightly different orientations, in a few cases it was 
not possible to completely exclude the intensity of these 
grains, which unfortunately impacted on the data qual- 
ity. 

Another problem that arises with the single-crystal 
data is twinning. Since the three-fold symmetry of the 
triangular planes is weakly distorted, three structural 
twins can exist rotated with respect to each other by 
60° in the be plane, while sharing the out-of-plane a-axis. 
Figure [6] shows all the possible nuclear reflections for all 
three twins in the (1, k, I) plane. Some reflections arise 
from one twin only e.g. (1, 1, 0); and can be used to deter- 
mine the ratio of the twins by simply taking the intensity 
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to the twinned reflections (each in the coordinate system 
of that twin). Using the atomic coordinates obtained 
from the E9 diffraction data at 2.1 K, this dataset was 
refined to obtain the scale factor and the volume ratios of 
the twins. The three twins were found to have non-equal 
weight (see Table , and their extracted ratios were es- 
sential to the magnetic structure refinement of the low 
temperature single crystal data which is described in the 
next section. 

TABLE II. Twin ratio of the single-crystal sample, deter- 
mined from single-crystal diffraction on E5 at 6 K. 



Twin 1 


Twin 2 (+60°) 


Twin 3 (+120°) 


64(2)% 


19(2)% 


17(3)% 



FIG. 5. Temperature dependence of the lattice parameters 
obtained from ID31 data. The lines are guides to the eye. 



C. Magnetic structure 



ratio of equivalent reflections of the three twins. Other 
reflections have overlapping contributions from all three 
twins which cannot be resolved on E5 e.g. (1, 0, 2) which 
is labeled (1, 2, 1) and (1, 2, —1) in the other two twins. 
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FIG. 6. The (l,k,l) reciprocal space plane where distortion is 
exaggerated for better visibility. The crosses give the position 
of the nuclear reflections while the dots show all possible mag- 
netic reflections. Each color correspond to a different twin. 
The dashed line indicates the border of magnetic Brillouin 
zone. 

A single intensity list was generated containing reflec- 
tions of all three twins and each reflection was labeled by 
the twin(s) to which it belonged. One index was assigned 
to the untwinned peaks while three indices were assigned 



Several different neutron diffraction techniques were 
necessary to determine unambiguously the magnetic 
structure of a-CaCr 2 04; these were powder diffraction, 
single-crystal diffraction and spherical neutron polarimc- 
try. We will discuss the results of each and show how 
together they provide a complete picture of the magnetic 
order. 

The neutron powder diffraction pattern from D20 mea- 
sured well below Tjy at 10 K is shown in Fig. [7] Several 
new peaks are observed which are attributed to the mag- 
netic order, the strongest is at d = 4.093 A where the 
broad bump was found in the 50 K diffraction pattern. 
The new peaks were indexed using the k-search program 
of the FullProf suite, the magnetic ordering wavevec- 
tor k was found to be close to the commensurate value 
k = (0, 1/3, 0) in agreement with Ref. [HI see Table [HI 
Only magnetic reflections where h is odd are observed, 
which implies antiferromagnetic ordering between suc- 
cessive planes along a (since the unit cell contains two 
Cr 3+ layers stacked along a). All possible magnetic re- 
flections in the (1, k, I) plane generated by k = (0, 1/3, 0) 
are shown in Fig. [6] However, only magnetic peaks 
where the three twins overlap have magnetic intensity 
e.g. (1,4/3,0); this makes the magnetic refinement chal- 
lenging. 



TABLE III. The magnetic peaks indexed assuming k = 
(0, 1/3, 0) are listed along with their observed and calculated 
d-spacings. 



dobs (A) 


index 


deal (A) 


4.093 


(l,l,0) + k 


4.084 


2.833 


(3,l,0)+k 


2.832 


2.154 


(l,3,0)-k 


2.154 


1.642 


(l,3,3)-k 


1.645 



The magnetic structures allowed for a-CaCr 2 04 can be 
derived from symmetry analysis, see Ref. 1201 A number 
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FIG. 7. Powder diffraction data measured on D20 at 10 K are 
displayed along with the Rietveld refinement assuming helical 
magnetic order in the ac plane, (a) The open red symbols are 
the data, the black curve is the fit, the upper and lower green 
bars indicate the nuclear and magnetic Bragg peak positions 
respectively, and the lower blue line gives the difference be- 
tween fit and data, (b) Detail of the (1, 1, 0)+k magnetic peak 
along with the 50 K data (open black symbols). The green 
and blue solid lines are the difference between the data and 
the fit with k fixed to (0,1/3,0) and refined to (0,0.332(3),0) 
respectively. 



of magnetic structures are possible including helical, cy- 
cloidal, ellipsoidal and sinusoidal with the spin moments 
pointing along different directions or rotating in differ- 
ent planes. Table [IV] gives the fitted magnetic mo- 
ments along the crystallographic axes. The refined phase 
between the Cr(l) and Cr(2) sites was close to — 2tt/3 , 
which was afterwards fixed to this value for all the refine- 
ments. The moment size on the two inequivalent Cr 3+ 
sites were refined independently, but found to have the 
same size within standard deviation. Therefore they were 
constrained to be equal for all subsequent refinements. 
As can be observed, the results of the refinement are am- 
biguous since several structures have similar i?-values. 
The powder data favors models where the spins point ei- 
ther in the ah or ac planes; however in these cases it is not 



possible to distinguish between helical or sinusoidal order 
because the phase difference between the two prefactors 
cannot be refined. Nevertheless the sinusoidal structures 
can probably be excluded indirectly as shown in Ref. |2T)1 
For example the refined moment lengths for the ac struc- 
ture are m x = 2.50(4)/ib and m z — 1.88(7)/xb- Assum- 
ing sinusoidally modulated order, the maximum moment 
length would exceed the spin only magnetic moment of 
Cr 3+ : yj m\ + m\ > gS 7 orbital momentum is quenched 
for Cr 3+ in an octahedral environment and g is the Lande 
g-factor (g rs 2 for spin-only magnetic moment) . On the 
other hand the maximum size of the moment in the he- 
lical structure would be the x component 2.50/is. 

The refined value of the k-vector was (0, 0.332(3), 0) 
from D20 at 10 K and (0, 333(5), 0) from E9 at 2.1 K. 
According to the refinement the k-vector is equal to 
(0,1/3,0) within standard deviation, the difference be- 
tween the diffraction patterns generated by these two 
k values is negligible (see Fig. [71(b)). Our results are 
also in agreement with the more precise incommensurate 
value of (0,0.3317(2), 0) published previously^ 1 . In the 
case of helical structures, the phase of —2tt/3 together 
with k = (0, 1/3, 0) means that the angle between all in- 
traplanar nearest neighbor spins is 120°. This structure 
is expected only for the ideal triangular lattice where 
all nearest neighbor exchange interactions are equal and 
Heisenberg. It is therefore surprising to find it on the 
distorted triangular lattice of a-CaCr 2 04. 

As a complementary study, we performed single-crystal 
diffraction at 6 K in the magnetic phase. The twin ratio 
and scale factor were fixed to the values obtained from 
the single-crystal nuclear peaks and the data were refined 
with the same set of possible magnetic structures as for 
the powder refinements. The i?p-factors and resulting 



magnetic components are listed in Table IV The single- 
crystal data is not able to distinguish between the helical 
and sinusoidal models either. The measurable |M(Q)| 2 , 
the absolute value square of the magnetic structure fac- 
tor, is identical for a helical structure and a sinusoidal 
structure with same moment components due to the or- 
thorhombic symmetry. In the case of ac sinusoidal struc- 
ture, two magnetic S-domains can exist in the sample, if 
the two domains have equal weight the resulting magni- 
tude of the structure factor is the same as in the case of 
the helical structure. For the plane in which the moments 
are aligned, the best fit gives the ac plane in agreement 
with the powder diffraction. 

The temperature dependence of two magnetic peaks 
(1,2/3,-1) and (1,4/3,0) were measured on the single- 
crystal sample in the temperature range 6 < T < 50 K, 
see Fig. [8j The intensity is proportional to the square of 
the magnetisation which is the magnetic order parame- 
ter and can be fitted with I(T) = C ■ ((T - T N )/T N ) 2/3 , 
where C is an overall constant and (3 is the critical ex- 
ponent. The data were fitted between 32 K and 48 K 
and gave a transition temperature of Tn = 42.6(1) in 
agreement with the heat capacity data. The exponent 
is j3 = 0.1815(9) for both peaks. Although the function 



TABLE IV. Refinement of the magnetic structure for spins pointing in various directions or planes for the two data sets (D20 
and E5). The column headings define the axes along which the magnetic moment is nonzero and m x , m y and m z are the 
moment sizes along these axes. All moment lengths are in units of zjb- The i?-factors are also included in order to compare the 
quality of the fits, Rb is determined from the integrated intensities: Rb = 2~Z I^obs — Jcalc|/X] I Jobs |- Helical and sinusoidal 
models give the same R-factors. 









ab 




ac 




be 


a 


b 


c 


D20 data (T = 10 K): 


m x 


2.62(5) 


2.18(1) 


2.50(4) 


2.24(2) 
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E5 data (T=6K): 


m x 


2.34(3) 


2.09(1) 


2.26(2) 
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2.99(2) 








m y 


1.78(4) 
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m x 
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3.14(2) 




Rf 


0.338 


0.341 


0.321 


0.322 


0.392 


0.392 


0.379 
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fits well, according to theory the magnetic phase transi- 
tion of a stacked triangular lattice is first order, and the 
value of (3 is not universal.^ Several noncollinear easy 
axis antiferromagnets have a sinusoidal phase just be- 
low Tn which precedes the helical phase, e.g. CuCrOj^ 
and /3-GaCr2 04.^ There is no sign of such a phase how- 
ever in a-CaCr 2 04, the temperature dependence of the 
two measured Bragg peaks is smooth and the ordering 
wavevector remains constant up to Tn (inset of Fig. ^ . 
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Temperature (K) 

FIG. 8. Temperature dependence of the (1,2/3,-1) and 
(1,4/3,0) magnetic reflections from the single-crystal sample, 
measured on E5. The black line is the fit described in the 
text. The inset shows the temperature dependence of k y re- 
fined from D20 data. 

Spherical neutron polarimetry was used to produce ad- 
ditional information to complement the previous neutron 
diffraction techniques. A number of different overlap- 
ping magnetic Bragg peaks were measured for all three 
twins. For example the strongest magnetic Bragg peak, 
(1,4/3,0) indexed in the coordinate system of twin 1, 
overlaps with (1, 2/3, 1) from twin 2 and (1, —2/3, 1) from 
twin 3. The (1, 4/3, 0) peak was also measured for twins 2 
and 3, where it is described as (1, 2/3, 1) and (1, —2/3, 1) 
respectively in the coordinate system of twin 1 (see Fig. 



[6]). The polarisation matrices for these three peaks la- 
beled using the coordinate system of twin 1 are displayed 
in the first three rows of Table [V] The observed data 
show strong depolarisation of the scattered beam when 
polarisation of the incoming beam was y or z due to the 
presence of structural twins. This makes it difficult to 
draw conclusions about the axis of the spin helix by sim- 
ply inspecting the data. Therefore, simulations were per- 
formed of the different models of magnetic ordering and 
compared to the data. The polarisation matrices were 
calculated using Blume-Maleev equations^ assuming the 
phase difference between Cr(l) and Cr(2) is — 27r/3 and 
the ordering wavevector is (0, 1/3, 0). 

Polarisation matrices of three possible spin structures 
(helical order in the ab, ac and be planes) were simulated 
for all three twins for the (1,4/3,0) magnetic reflection 
in the coordinate system of twin 1, see Table [yi] For the 
helical models the issue of chirality also needs to be ad- 
dressed. In these cases the centrosymmetry of the param- 
agnetic space group is broken in the magnetically ordered 
phase, allowing the appearance of two chiral domains for 
each twin with opposite vector chirality. The polarisa- 
tion matrix of the chiral domains differs only in the V yx 
and "P zx terms, which change sign with opposite chirality, 
see Table IVT1 Because the measured values for these two 
elements are rather small (less than 0.136(9)), it is as- 
sumed that the population of the two chiral domains are 
equal for all three twins. Fitting the weight of the chiral 
domains also don't improve the data significantly. In this 
case the measured V yyi , V zx would be zero because they 
cancel each other out. This, together with the presence of 
structural twins, explains the strong depolarisation when 
the incident beam has y and z polarisation. In order to 
make a comparison with the measured (1,4/3,0) reflec- 
tion, the contributions of the three twins weighted by 
their volume ratios were added together and multiplied 
by the beam polarisation and are listed in the last three 
columns of table |Vl[ Comparison of the simulated matri- 
ces to the experimental matrix in the first three rows of 
Table |v| under (1, 4/3, 0), clearly show that the ac helical 
structure gives the best description of the data. 

To precisely determine the orientation of the spin ro- 
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TABLE V. Polarisation matrices measured using TASP at 1.5 K for the three strongest magnetic reflections (rows 1-3) and 
results of the refinement which yields a helical structure with spins lying almost entirely within the ac plane (rows 4-6). 
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0.372(8) 
0.082(9) 


0.033(10) 
0.063(9) 
-0.369(9) 


-0.903(4) 
-0.001(9) 
-0.014(9) 


-0.036(10) 

0.400(9) 
-0.092(9) 


0.004(10) 

0.110(9) 

0.402(9) 


px 


-0.887 


0.000 


0.000 


-0.887 


0.000 


0.000 


-0.887 


0.000 


0.000 


Pf 


0.000 


0.098 


-0.022 
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Pf 
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TABLE VI. Simulations of the polarisation matrices of the (1,4/3,0) reflection for helical order in the ab, ac and be planes. 
The ± sign is for the two opposite vector chiral domains. The contributions of all three structural twins are listed. Columns 
10-12 give the combined contribution weighted by the ratio of the twin volumes and assuming equal chiral domain ratios, the 
last three columns are multiplied by the finite instrumental polarisation. 



Twin 1 Twin 2 (+60°) Twin 3 (+120°) Together Pxx corrected 
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Pf Pf Pf 
-0.883 0.000 0.000 
0.000 0.496 0.005 
0.000 0.005 -0.496gl 


px 

ac Pf 
Pf 
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FIG. 9. Magnetic structure of Q-CaCi'204, determined from 
spherical neutron polarimetry. Cr(l) and Cr(2) are repre- 
sented by the blue and red spheres respectively. 



tation plane, the polarisation matrices of several other 
reflections were also considered: (3,4/3,0), (5,4/3,0), 
(1, 8/3, 0) for all three twins. A computer code was writ- 
ten to fit these matrices to both sinusoidal structures with 
arbitrary spin direction and helical structures with arbi- 



trary spin rotation plane. Spherical coordinates (9, if) 
were used to parametrize the spin direction in the case of 
sinusoidal structures, or the vector n normal to the ro- 
tation plane in the case of helical structures (where the 
c axis corresponds to 9 = 0°). The equally populated S- 
domains were generated from the magnetic structure by 
applying one of the mirror planes (m yz and m xz ) or the 
glide plane (n xy (l/2, 1/2, 0)) and the contribution from 
all domains were averaged. The other fitted parameters 
were the beam polarisation (V xx ) and the ellipticity of the 
spin helix (ratio of the moment length along the two main 
axis of the ellipsoid, where /i = m v /m u , u = n x (c x n) 
and v = ux n). Fitting was performed by minimizing the 

function R%, = X>|"Pobs ~ ^calc| 2 / E ^obs- where 
w = 1/(T 2 and a is the counting error and the summa- 
tion was over all elements of all measured polarisation 
matrices. 

By fitting this general model to all the measured po- 
larisation matrices the magnetic structure was found to 
be helical with spins rotating in the ac plane; the result- 
ing polarisation matrices are listed in the last three rows 
of Table |V| The extracted parameters are: "P xx = 0.887, 
9 = 87°, tp = 91°, m a /m c = 0.99. With these parame- 
ters excellent agreement with the data could be achieved 
(R w = 0.074). To check how sensitive the refinement is 
to the spin-plane orientation, the values of 9 and (f were 
fixed to other directions. All other possible structures 
(ab and be helical, a, b and c sinusoidal) give a value of 
R w that is more than four times larger. According to 
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the result the spin components along a and c are close 
to equal in agreement with Ref. [20l whereas the single 
crystal and powder refinements show larger a component 
by a factor of 1.2. 

Using spherical neutron polarimetry it has been pos- 
sible to prove unambiguously that a-CaCr 2 04 has a 
~ 120° helical magnetic structure with spin moments 
rotating in the ac plane perpendicular to the ordering 
wavevector k = (0, ~ 1/3,0). The magnetic structure is 
illustrated in Fig. [9] 



IV. DISCUSSION 

The 120° helical magnetic structure of a-CaCr 2 04 is 
exactly the structure expected for an ideal triangular 
antiferromagnet with Heisenberg nearest neighbor inter- 
actions that are all equal. This is a surprising result 
given that the triangular plane is in fact distorted with 
two inequivalent Cr 3+ ions and four independent near- 
est neighbor exchange interactions. Other distorted tri- 
angular antiferromagnets that have been investigated so 
far all show departures from ideal triangular behaviour 
with ordering wavevectors that deviate from (0,1/3,0) 
e.g. Cs 2 CuCl 4 (where k = (0,0.472,0)), NaMn0 2 (where 
k = (1/2,1/2,0)) and CuCr0 2 (where k = (0,0.329,0)). 
In this section we solve this apparent contradiction by 
investigating how the ordering wavevector and relative 
orientation of the spin moments in a-CaCr 2 04 vary as a 
function of the exchange interactions. 

The possible interactions between the magnetic mo- 
ments can be deduced from the crystal structure. Figure 
[j] shows the location of the Cr 3+ ions on the triangular 
plane. There are four different nearest neighbour Cr 3+ - 
Cr 3+ distances, d c hi, d c h2, d zz i and d zz2 , ranging from 
2.882 to 2.932 A at 2.1 K which are identified by the 
different colored lines. The Cr 3+ ions are surrounded by 
oxygen octahedra which lower the energy of the £ 2g or- 
bitals and since there are three electrons in the 3d shell, 
each of the t 2g orbitals is occupied by one single elec- 
tron. Neighboring CrOg octahedra are edge-sharing and 
this along with the short distances between the magnetic 
ions implies that the exchange interactions are dominated 
by direct overlap of singly occupied i 2g orbitals and are 
therefore antiferromagnetic!2SI Since direct exchange in- 
teractions are highly sensitive to ionic separation we ex- 
pect four different exchange constants, J c hi, J c h2, Jzzi 
and J zz2 , corresponding to the four nearest-neighbor dis- 
tances as labeled in Fig. [T] These interactions are ex- 
pected to be Heisenberg with no anisotropy because the 
occupied t 2g multiplet has quenched orbital angular mo- 
mentum. Although there are two triangular layers per 
unit cell, they are in fact the same due to the the reflec- 
tion plane m(l/4, y, z) of the Pmmn space group. The 
interlayer interactions are also expected to be antifcr- 
romagnetic but much weaker and to occur via superex- 
change through intermediate oxygen ions. Furthermore 
these interactions are unfrustrated and do not influence 



the ordering k wavevector. Next nearest neighbor intra- 
plane interactions may be present and would influence 
k, as reported for CuCr0 2 , 17 another member of the 
delafossite arisotype which has similar Cr 3+ -Cr 3+ dis- 
tances. 




FIG. 10. Magnetic phase diagrams as a function of Ai, A2, 
and Jnnn with A3 fixed to zero for all graphs. Inside the blue 
region k = (0, 1/3, 0) with 120° spin structure ((p la = 0, (p lb = 
240°, ip 2a = 120° and <p 2b = 0°). Elsewhere k = (0, 1, k z ) and 
the spin angles vary. 



The magnetic Hamiltonian of a-CaCr 2 04 can be ex- 
pressed as the sum of a Hamiltonian containing nearest- 
neighbor interactions and a Hamiltonian containing next- 
nearest-neighbor interactions, while the much weaker in- 
terlayer interactions have been ignored 
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The nearest neighbor Hamiltonian has the form 



into the results as will be explained later 
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where S nm is the spin of Cr 3+ site (i<E {la, lo, 2a, 2o}), 
and n and m are the indices of the unit cell along the 
b and c axes respectively. The next-nearest-neighbor 
Hamiltonian has the form 
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where a single J nn n exchange constant is assumed for all 
next-nearest-neighbor interactions. 

The Hamiltonian can be used to explore the magnetic 
structure of a-CaCr 2 04 as a function of the exchange in- 
teractions. In order to do this all the spins are treated as 
classical vectors of the same length and are restricted to 
be coplanar. The aim is to understand how the in-plane 
k-vector (both k y and k z components) and the directions 
of the spin moments (p nm (for the ith spin of the (n, m)-th 
unit cell) depend on the magnetic interactions. </? nm can 
be written in terms of if 1 , the angle of spin S l in the first 
unit cell, and the ordering vector: (^ nm = (p l + k • r nm , 
where r nm is the position of the unit cell. The Hamilto- 
nian can be re-expressed in terms of k and (p 1 by rewriting 
the dot product of the spins, for example 
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[{Jchl + Jch2) — (Jzzl + Jzt.2)]- (5) 



■/mean gives the overall energy scale and does not affect 
the zero temperature ground state structure, thus to- 
gether with J nnn only four important parameters remain. 
Finally, to determine the ground state values of k, ip lb , 
tp 2a and (f 2b for a given set of exchange interactions Ai, 
A3 , A3 and Jnnn the energy of the Hamiltonian was min- 
imized by performing a constrained linear optimization 
using a simplex minimisation algorithm. 



Figure 10 shows the in-plane k-vector as a function 
of Ai, A 2 and J nn n with A 3 fixed to zero. The reason 
for the linear transformation of the exchange constants 
is clear from these plots. They show that the magnetic 
structure is commensurate (k = (0, 1/3,0)), with a 120° 
ground state (tp la = 0, tp lb = 240°, tp 2a = 120° and 
ip 2b = 0° over a large volume of Ai, A 2 and J nnn param- 
eter space. Although the crystal structure is distorted, 
the nature of the distortion is such that when A3 is con- 
strained to zero, commensurate 120° magnetic order is 
stabilized unless the magnitude of Ai, A 2 and J nnn are 
much greater than 0. 
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ip la can be set to zero due to the 0(3) invariance of the 
Heisenberg Hamiltonian. 

A linear transformation was performed on the ex- 
change constants, which was found to give more insight 



FIG. 11. Magnetic phase diagrams showing the dependence 
of k on Ai, A2 and A3; (a) displays k y and (b) k z in units 
of 2n. For A3 = 0, k = (0,1/3,0) and the spins form a 
120° spin structure (<p la = 0, ip lb = 240°, ip 2a = 120° and 
ip 2b — 0°) within the red border. For A3 7^ 0, k becomes 
incommensurate . 



If A 3 is varied as well the k-vector changes contin- 



uously as shown in Fig. 11 Within the red borders on 
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the plots, all four ip l angles are determined only by the It- 
vector, according to the expression ip 1 = (k+ (0, 1, 0)) • p 1 
where p 1 is the position of S 1 in the first unit cell. Out- 
side these regions the ordering wavevector locks into a 
new set of value (k = (0, 1, k z )) and the relative phases 
inside the unit cell are independent of k. Figure |l2"| 
shows how k varies with A 3 , (all other parameters are 
fixed to zero). Even a slight deviation of A3 from zero 
drives k incommensurate. The slope of k y at the point 
A3 = is —0.74, this can be used to estimate the value 
of A 3 for a-CaCr 2 04. From powder diffraction data 
the upper limit for the difference between k y and the 
commensurate 1/3 value is d(k y ) < 0.004. This means, 
|A 3 | < |l/(-0.74)|d(fc y ) - 0.005. So the only condition 
for the observed magnetic structure is that |A 3 | < 0.005. 
Comparison of the different Cr 3+ -Cr 3+ nearest neighbor 
distances reveal that the average of d zz i and d ZZ 2 (2.907 
A) is almost the same as the average of o? c hi and d c h2 
(2.912 A). This explains why the direct exchange inter- 
actions appear to obey (J zz i + J zz2 ) ~ (J c hi + ^cha) giving 
|Aa|«0. 




-1.0 -0.5 0.0 0.5 1.0 



A 3 

FIG. 12. Dependence of the y component of k on A3 
(Ai = A2 = Jnnn = 0), right scale shows the phase differ- 
ence between S la and S 2a . At A 3 = the tangent of the 
slope (blue line) is -0.74. 

Although a-CaCr2 04 is distorted from ideal triangu- 
lar symmetry with two inequivalent Cr 3+ ions and four 
different exchange interactions, the nature of the distor- 
tion is such that the average of the exchange interactions 
along any direction are approximately equal. On larger 
length scales the Hamiltonian is in fact spatially isotropic 
providing an explanation for the observed ~ 120° helical 
magnetic ordering. 



V. CONCLUSION 



In this paper we presented a detailed investigation of 
the magnetic order in the distorted triangular antiferro- 
magnet a-CaCr2 04. The first single-crystal growth is 
reported and bulk properties measurements along with 
powder and single-crystal X-ray and neutron diffraction 
were performed. The ordering wavevector was found to 
be close to commensurate k = (0, 1/3,0) and we unam- 
biguously determined that the magnetic ordering is he- 
lical with the spin moments forming a ~ 120° structure 
between nearest neighbors in the ac plane. It is impor- 
tant to note that a small deviation from the commen- 
surate wavevector is probable due to the orthorhombic 
crystal symmetry. The apparent contradiction between 
the observed ideal magnetic structure and the distorted 
crystal structure was resolved by exploring the magnetic 
phase diagram. The distortions in the triangular plane 
are only significant on individual nearest neighbor dis- 
tances, while on the length scale of a unit cell the aver- 
age exchange interactions along any triangular axis are 
approximately equal. 
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